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One of the central problems of the liquid-glass transition is whether there is a structural signature 
that can qualitatively distinguish different dynamic behaviors at different degrees of supercooling. 

Here, we propose a novel structural characterization based on the spatial correlation of local density 
and we show the locally dense-packed structural environment has a direct link with the slow dynamics 
as well as dynamic heterogeneity in glass-formers. We find that particles with large local density 
relax slowly and the size of cluster formed by the dense-packed particles increases with decreasing the 
temperature. Moreover, the extracted static length scale shows clear correlation with the relaxation 
time at different degrees of supercooling. This suggests that the temporarily but continuously 
formed locally dense-packed structural environment may be the structural origin of slow dynamics 
and dynamic heterogeneity of the glass-forming liquids. 

PACS numbers: 61.20.Ja, 61.20.Lc, 64.70.P- 


Upon cooling, glass forming liquids display markedly 
slow and heterogeneous dynamics but unperceivable 
change in its static pair correlations 0, i- Although 
intense attention has been paid to explore the structural 
origin of slow dynamics and dynamical heterogeneity, a 
more convincing view is still lack and the relationship 
between structure and dynamics seems to be quite com¬ 
plex ii- It is still unclear whether the glass transition 
is due to an underlying thermodynamic phase transition 
or purely dynamic process in nature. The related mode¬ 
coupling theory, which predicts the dynamics solely based 
on the static pair correlations, is found to be failed to ex¬ 
plain the nonperturbative effect of attractive forces in su¬ 
percooled liquids, suggesting the inefficiency of pair cor¬ 
relation when describing the structural change of glass¬ 
forming liquid [m. 

To unveil static quantities, several different ideas have 
been proposed. The locally preferred structural order, 
accompanied with little density change, is claimed to 
be the cause of the slow dynamics, as growing in its 
length scales with increasing the degree of supercooling 
in several systems However, this method shows 

clearly system dependences [12|, and is hard to be ex¬ 
tended to different glass-formers (e.g., the high dimen¬ 
sional liquids [13 and the strong frustrated glass form¬ 
ers M)- Recently, a more general method investigating 
the point-to-set (PTS) correlations in the amorphous sys¬ 
tem emerged and attracted lots of attention It 

is worth noting that the obtained modest increase of the 
static PTS length scale (^pts) shows clear decoupling 
with the dramatic increase of the dynamical length char¬ 
acterizing spatial heterogeneities of the glass-formers [ 13 - 
[2lj , which arouses heated debates on the relationship be¬ 
tween structure and dynamics of glass-formers. Clearly, 
exploring more general structural signatures of slow dy¬ 
namics is crucial to understand the structural origin of 
the slow dynamics and dynamic heterogeneity for glass¬ 
forming liquids. 


For glass formers, one of the most pronounced features 
is that particles are trapped in the transient cages 22- 


25l| . which leads to the rigid nature of amorphous ma¬ 


terials. This brings an intriguing question that whether 
this localization dynamics corresponds to the change of 
the local structural environment. A perfect candidate, i. 
e., the local density (p/), which can be determined from 
the Voronoi cell volume (v), can provide such local infor¬ 
mation for each particle. Previous investigations found 
weak correlation with dynamics based on per-particle lo¬ 
cal geometry mill. Up to now, neither indications of 
the spatial density correlation nor the link between the 
local structure characterized by pi and glassy dynamics 
was reported. 

In this work, instead of concentrating on the first 
coordination shell local geometry, we study the spatial 
correlation of the relatively dense-packed structural en¬ 
vironment and explore its links with the slow dynam¬ 
ics as well as dynamic heterogeneity in glass-forming 
liquids by means of molecular dynamics simulation. 
We achieve positive results of the above relationship 
by studying three dimensional (3D) Kob-Andersen bi¬ 
nary Lennard-Jones mixture (LJ) |28| a nd its Weeks- 
Chandler-Andersen (WCA) truncation [29| over the den¬ 
sity range from p = 1.2 to p = 1.9 (see system details 
in the Supplemental Material Hi). These two systems 
have different dynamics at low densities in the super¬ 
cooled regime while this difference shrinks as increasing 
p and almost disappears at p = 1.9 0 0 113 • On the 
other hand, in the whole studied density range, the two 
systems have nearly identical two body static correlations 
at each temperature. Thus these two systems provide a 
perfect testbed to study whether the local density is sen¬ 
sitive to structural differences and whether there exists 
a direct link between structure and dynamics for glass- 
formers. We demonstrate that LJ system has a larger 
size of clusters (< Nc >) formed by the relatively dense- 
packed particles than WCA system at lower density. This 
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FIG. 1: (color online). Scaled distribution of Voronoi cell 
volumes, avP{v) as a function of the normalized Voronoi vol¬ 
ume, ((u— < V >)/(T„) at 30 different state points of variable 
T and p (ranging from high T liquid to low T supercooled 
states for both LJ and WCA systems). All the results are for 
A particles. 


corresponds to the slower dynamics of the LJ system. 
Moreover, the difference of the value of < A"c > for these 
two systems shrinks synchronously with the dynamics as 
increasing p. The extracted static length scale is found to 
correlate well with the degree of supercooling and also the 
relaxation time in the studied temperature and density 
range, suggesting that the locally dense-packed environ¬ 
ment might be the main reason of slow dynamics and 
dynamical heterogeneity in glass-formers. 


To find which particles are densely packed in the sys- 
tem, we first implement the Voronoi constructions 


and calculate the Voronoi cell volume v for each particle. 
Follow the protocol of Ref. [s^, we investigate the sta¬ 
tistical properties of the Voronoi cells by calculating the 
scaled distribution of v (see Fig. 1). Here, v is shifted by 
the average Voronoi volume < v > and then normalized 
by the standard deviation (t„ (cr^ =< > — < v >^). 

As expected, the distributions {P{v)) only show little dif¬ 
ference for various T and p and are not sensitive to the 
system difference, consistent with Refs, Such 

kind of scaling collapse behavior is expected to be uni¬ 
versal in dense liquids with relatively harder interactions, 
and deviations have been found in soft-potential colloidal 
suspensions [dlj. Our results indicate similar intrinsic 
amorphous structure from the statistical view at differ¬ 
ent T and p of the LJ and WCA glass formers. Although 
neither the static pair correlations nor the statistical dis¬ 
tribution of the Voronoi volumes accounts for the sig¬ 
nificant dynamic differences of these two systems, other 
subtle characterizations such as the higher order struc¬ 
tural correlations 11, 4^, the PTS correlation length 181 
and properties of zero-temperature glasses can also 
capture the structural differences. 

The distributions {P{v)) show little dependence on 
T and p for both LJ and WCA systems, however, the 
growth of the static correlations of the relatively dense- 


packed particles exhibit marvelous connection with slow 
dynamics upon cooling. First, we emphasize that the lo¬ 
cal density is not correlated with local ordering, as has 
been proved in Refs. [ 1 , We define 13% of parti¬ 

cles with the largest local density pi {pi = 1/v) as the rel¬ 
atively dense-packed particles (The qualitative results are 
not sensitive to the choice of the threshold value 13%). 
We also use other criterions to define the densely packed 
particles (e. g., particles with Voronoi volume satisfying 
{v— < V >)/(jy < —1.2 are defined as the densely packed 
ones) and we can obtain the very similar results (see the 
Supplemental Material [s^ for details). For type A and 
type B particles, the range of pi is distinct, thus we sort 
type A and type B particles separately and choose 13% 
of type A particles and 13% of type B particles with the 
largest local density respectively as the densely packed 
particles. Two densely packed particles are considered 
as belonging to the same cluster if they are neighbors 
defined by the Voronoi method. The number of densely 
packed particles that are belonging to a cluster is defined 
as the single cluster size Nc- The number averaged clus¬ 
ter size < Nc > is defined as 


< Nc >= 


J:Z=i NcPiNc) 


( 1 ) 


where P{Nc) is the probability of finding a single cluster 
size Nc, and = 1- The value of < TVc > can 

reflect the degree of spatial correlations between dense- 
packed particles. 

In Fig. 2, we show the temperature dependence of the 
relaxation time Tq, and that of the number averaged clus¬ 
ter size < Nc > at p = 1.2 and p = 1.9 for both LJ 
and WCA systems. Ta is defined as Fs{qp,t = To) = 1/e 
(where Fs{qp, t) is the intermediate scattering function of 
A particles as shown in the Supplemental Material [23). 
It is clear that both Tq, and < Nc > increase as the tem¬ 
perature is lowered, implying that < Nc > can reveal un¬ 
ambiguous structure difference between the supercooled 
liquid and normal liquid as an effective static signature. 
At p = 1.2, LJ system has larger Tq, and thus its dy¬ 
namics is slower than the WCA system at the same T in 
the supercooled states (see Fig. 2(a)). Correspondingly, 
< Nc > is larger for LJ system (see Fig. 2(c)). On the 
other hand, the < Nc > difference of the two systems 
is small but far from negligible at high T, although Tq, 
of the two systems is quite close. This is quite reminis¬ 
cent of the predictions from the Arrhenius law that the 
LJ and WCA system have distinct dynamics in the high 
T regime at p = 1.2, since the active energy of the two 
systems is different [y]. At p = 1.9, both and < Nc > 
are very close for the two systems, implying both dynam¬ 
ics and its structural signature are similar (see Fig. 2(b) 
and 2(d)). Our results show direct evidence of the cor¬ 
relation between static property characterized by the lo¬ 
cally dense-packed structural environment and dynamic 
behavior at different degrees of supercooling in a large 
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FIG. 2: (color online). T dependence of the relaxation time 
Ta of the LJ and WCA systems for p = 1.2 (a) and p = 1.9 
(b), and T dependence of the number averaged cluster size 
< A^c > of the LJ and WCA systems for p = 1.2 (c) and 
p = 1.9 (d). The black squares and red circles are for LJ and 
WCA systems, respectively. < Nc > is obtained by averaging 
over 1000 configurations in each run, and the error bars in (c) 
and (d) are obtained by averaging over 8 independent runs. 


density range. It should be noted that, different from the 
static PTS correlation length the present < > 

data undergoes no fitting or further processing, and this 
unambiguously avoids some sorts of artifacts in the in¬ 
terpretation of the data. This local density based spatial 
correlations can be an overwhelming tool to characterize 
the static nature of glass-formers in some systems such 
as the high dimensional liquids [l^ or strong frustrated 
glass-formers 1J|, in which the bond-orientational order 
is not obvious. 

We define the static correlation length of locally dense 
structural environment as ^ Nc> based on our 

three dimensional glass model. Similar definition in two 
dimensional glass-forming liquids has been used in defin¬ 
ing the length scale of the hexatic structural order [i€3- 
Next, we study the growth of ^ as increasing the degree 
of supercooling and its link to slow dynamics for all the 
density range we studied in both LJ and WCA systems. 
Fig. 3(a) illustrates the behavior of Tq as a function of 
the static length scale It is clear that Ta grows mono- 
tonically with ^ in the whole temperature range for all the 
systems, which is consistent with the behavior of other 
reported static length scales [S-T^, 16-21| as increasing 
the degree of supercooling. Moreover, the data in Fig. 
3(a) shows that the difference between LJ and WCA sys¬ 
tems is pronounced at p = 1.2, but shrinks as increasing 
p, indicating that attractive tails have a remarkable im¬ 
pact on the system at relatively low density, similar to 
the related reports [11,1^. Inspired by the temperature de¬ 
pendence of the length scale of the structural order Q, 
we expect that ^ behaves as 


^ = m{T-To)Y 


( 2 ) 



(l/Dlog(TyT„))^ 


FIG. 3: (color online), (a) The a relaxation time Tq vs the 
static length scale ^ for all the density range studied for both 
the LJ and WCA systems, (b) The rescaled length scale 
^/Co versus [ro/(r — To)]^. The black solid line is the re¬ 
lation of = [To/(T — Tb)]^. (c) Relation between 
and [1/Dlog{ra/To)]^■ The black solid line is the relation of 
^/Co = [1/Dlog{Ta/ra)]^■ Legends in (a), (c) are the same as 
those in (b). 


where and A are the fitting parameters. Both the 
value of them reflect the intrinsic property of the sys¬ 
tem and are quantitatively impacted by the definition of 
dense-packed particles. Detailed discussions o n fn and A 
can be found in the Supplemental Material [30|. Tq is 
the critical glass transition point fitted by Vogel-Fulcher- 
Tamman (VFT) relation 

Ta = Toexp[DTo/{T - To)], (3) 

where tq, D and Tq are the fitting parameters. D is the 
fragility parameter, and smaller value of D corresponds 
to a more fragile glass former. All our simulation data 
can be well fitted by VFT equation (see Fig. S2 of the 
Supplemental Material 0). Fig. 3(b) plots the relation 
between static length scale and temperature of the LJ 
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FIG. 4: (color online). Time dependence of the self part of 
dense-packed particles overlap and that of loose-packed 
particles overlap q^t) at T = 0.44 and p = 1.2 for both LJ 
and WCA systems. 


and WCA systems in the whole density range we stud¬ 
ied. Interestingly, all curves collapse when we plot 
against [To/(T—T q)]^, implying a direct link of the static 
length scale ^ and the degree of supercooling. Combining 
Eq. (2) and Eq. (3), we obtain the scaling relation of the 
static length scale ^ and the relaxation time Tq, 

i = U^/Dlog{TjTq)]^, (4) 

The results of this relation at different densities in the LJ 
and WCA systems are shown in Fig. 3(c). Clearly, the 
^ — Tq, relation supports the scaling argument of Eq. (4). 
As a result, our findings provide useful clues to the un¬ 
derlying causes of the slow dynamics in the supercooled 
regime. The static length scale of the densely structural 
environment plays a key role in determining the dynam¬ 
ics of the glass-forming liquids. Moreover, we test our 
results in other two glass models and we can also find 
similar collapses in the relations of both ^ — T and ^ — Tq 
( see the last section of the Supplemental Material [^). 
To our knowledge, this is the first finding of the link be¬ 
tween density based spatial correlations and dynamics. 

Another crucial question is the relation between the 
locally structural environment and the dynamic hetero¬ 
geneity. To characterize the dynamics of the dense- 
packed and loose-packed particles, the self part of dense- 
packed particles overlap q‘^{t) and loose-packed particles 
overlap q\t) are defined as follows. We partition the 
simulation box into small cubic boxes with side length 
Ig = 0.36, so that there is no more than one particle 
in each small cubic box. Then we define a binary digit 
nfit) to specify whether the same dense-packed particle 
or loose-packed particle occupies the cell i at times 0 and 
t. 


^ Ez < <W<(o) > 
Ez < <(o) > 


( 5 ) 


where the sum runs over all the small cubic boxes in the 
system. The superscript a = dis for dense-packed parti¬ 
cles and a = I is for loose-packed particles, respectively. 
Here, we also define 13% of particles with the largest 
local density as the dense-packed particles and 13% of 
particles with the smallest local density as loose-packed 
particles. 

In Fig. 4, we show the time dependence of q‘^{t) and 
q\t) at T = 0.44 and p = 1.2 for both LJ and WCA sys¬ 
tems. Similar to the self-intermediate scattering function 
as shown in Fig. SI [s^, the overlaps show clear two-step 
decay for both LJ and WCA systems, indicating that 
both two systems are in the supercooled regime. More¬ 
over, the LJ system has a higher shoulder for both g'^(t) 
and g^(t), which indicates its dynamics is slower than the 
corresponding WCA system. This is also consistent with 
our previous results. More importantly, the dense-packed 
particles relax more slowly than the loose-packed parti¬ 
cles by comparing q‘^{t) and q^{t) of the LJ system and 
WCA system. This phenomenon is more pronounced at 
lower temperatures (data of p = 1.2 for several tempera¬ 
tures for both LJ and WCA systems is shown in Fig. S3 
of the Supplemental Material [1^). This is reasonable 
since dense-packed particles may have less free volume 
to relax than the loose-packed particles. This also sug¬ 
gests that there exists direct link between local density 
and dynamic heterogeneity for glass-forming liquids. The 
relatively dense-packed particles relax slowly and domi¬ 
nant the slow relaxation dynamics. 

Our results indicate that the transiently localized 
cages, which trap the dynamics of particles, may have 
a length scale much larger than the size of fist coordina¬ 
tion shell, and growing in its size as further supercooled. 
The relatively dense-packed structural environment may 
be the structural presentation of low potential energy do¬ 
mains, since it has been found that particles with low po¬ 
tential energy move slowly and clusters formed by them 
correlate with the dynamic heterogeneity [i^. These 
findings support that the liquid-glass transition is not 
solely a dynamic behavior but also accompanied with 
subtle changes of structures. 

In Summary, we present evidences of the intrinsic link 
between the microscopic local density and the slow dy¬ 
namics as well as dynamic heterogeneity of the glass¬ 
forming liquids. We find, at a large density range, that 
the difference of the dynamics of the LJ and WCA sys¬ 
tems can be completely captured by the cluster size 
characterized by the relatively dense-packed particles, al¬ 
though the normally used static pair correlations and the 
statistical property of the Voronoi volume remain blind 
to the dynamics. The extracted static length scale of 
the dense-packed particles directly corresponds to the de¬ 
gree of supercooling and also the relaxation time of the 
glass-formers. The temporarily but continuously formed 
locally dense-packed structural environment may be the 
origin of the slow dynamics and dynamic heterogeneity in 
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the glass-forming liquids. Our findings open a new way 
to study the central problem of glass transition and may 
also stimulate more experimental work on the structural 
nature of slow dynamics and dynamic heterogeneity. 
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The system information 

The present studied systems consist of = 3074 — 
4865 (depending on p) with 80% A and 20% B par¬ 
ticles. The potential between particles a and (3 is 
Vap{r) = ‘3eap[{cTap/rapY‘^ “ {papjTa.pf + C^p], when 
r^p is smaller than the potential cutoff and zero 
otherwise. Here, a, f) G {A,B}. The interaction param¬ 
eters are given by oabI<^aa = 0.8, <jbbI<^aa = 0.88, 
^abI^aa = 1-5, and cbbI^aa = 0.5. The poten¬ 
tial is truncated and shifted at i.e., 2^/^aap for 
WCA system and 2.5 (Tq^ for LJ system. Cap guarantees 
Vap{r‘f,^p) = 0. Length, energy and time are recorded in 
units of (7aa, ^aa and y/ma\j^j€AA, respectively. The 
simulations are performed in the NVT ensemble with 
periodic boundary conditions. The Newton’s equations 
of motion were integrated with the velocity form of the 
Verlet algorithm and the temperature T was maintained 
by the Nose-Hoover thermostat UM- In order to ob¬ 
tain reliable results and improve the statistics, for each 



FIG. S 1: (a) The radial distribution functions (results have 
been shifted for clarity) and (b) the self-intermediate scat¬ 
tering functions at different temperatures for p = 1.2 (left) 
and p — 1.9 (right) for both LJ (black solid lines) and WCA 
(pink dashed lines) systems. The horizontal dotted line in 
(b) marks Fs{qp,t) — 1/e at which the relaxation time Ta is 
determined. All the results are for A particles. 



FIG. S 2: r dependence of the relaxation time Ta- Data are 
shown in the VFT form log(Ta/To) ~ DTo/{T — To), for all 
the temperatures at different densities of the LJ and WCA 
systems. The black solid line is the relation of log(Ta/To) = 
DTq/{T-To). 

state point, we average the results over eight independent 
runs. All the simulations are performed with in-house 
GALAMOST software package Q. 

We measure the intermediate scattering functions of 
A particles, F,{qp,t) = ^ < Ej^iexp{iQp • [rj{t) - 
'^j(0)]} >: where Vj is the position vector of particle j, 
< ■ • • > indicates the thermal average, i = '/-A and the 
wave number qp = 7.3 corresponds to the first peak of the 
static structure factor. The similarity of the static pair 
correlations and differences of the relaxation dynamics 
at different degrees of supercooling for p = 1.2 and p = 
1.9 are shown in Fig. SI. Clearly, the pair correlations 
cannot account for the differences in dynamics for glass¬ 
forming liquids. 

We show in Fig. S2 the VFT (Eq. (3) in the paper) fit¬ 
ting results of temperature dependence of the relaxation 
time in the whole density range we studied for both LJ 
and WCA systems. Clearly, the T dependence of Tq re¬ 
markably agree with the VFT equation for all the state 
points investigated. 

The comparison of q‘^{t) and g*(t) at different temper¬ 
atures for p = 1.2 in both the LJ system and the WCA 
system can be seen in Fig. S3. The curves show one- 
step decay at higher temperatures and two-step decay 
at lower temperatures, which is akin to the T depen¬ 
dence of the intermediate scattering functions. The q'^{t) 
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FIG. S 3: Time dependence of the self part of densely packed 
particles overlap (black solid lines) and that of loosely 

packed particles overlap q\t) (red dashed lines) for several 
temperatures at p = 1.2 for both the LJ system (a) and 
the WCA system (b). Temperatures for the LJ system are 
3.0, 1.0,0.6,0.5,0.48 and 0.44 from left to right in (a). Tem¬ 
peratures for the WCA system are 3.0,1.0, 0.5, 0.4, 0.3 and 
0.28 from left to right in (b). 


curves always lie up the curves of q^{t) before decaying 
to 0 even at high temperatures. This shows that the 
locally dense-packed particles always relax more slowly 
than the loosely packed particles, and this phenomenon is 
more pronounced in the supercooled regime. Results for 
the high p data show similar behavior (data not shown). 
Thus, it can be concluded that the dynamic heterogene¬ 
ity has a clear structural origin, i.e., the heterogeneity of 
local density. 



FIG. S 4: T dependence of the number averaged cluster size 
< Nc. > of the LJ and WCA systems for p = 1.2 (a) and 
p = 1.9 (b). < A^c > is obtained by averaging over 1000 

configurations in each run, and the vertical error bars are 
obtained by averaging over 8 independent runs. Here the 
definition of densely packed particles is different from that in 
the main text (see SM text). 


inition, we calculate the averaged cluster size < Nc > 
and show in Fig. S4 the < Nc > data of p = 1.2 and 
p = 1.9 for both the LJ and WCA systems. Compared 
with results shown in Fig. 2 in the main text, the values 
of < Nc > are smaller since the number of densely packed 
particles approximately accounts for 10% of the total par¬ 
ticles and this percentage is less than 13% defined in the 
main text. Despite of quantitative differences, the qual¬ 
itative results are very similar for both these two defini¬ 
tions of dense-packed particles. At p = 1.2, < Nc > is 
larger for LJ system and curves show big difference from 
WCA system. However, at p = 1.9 the two curves are 
very close, which implies that < Nc > has a clear con¬ 
nection with the dynamics, as has been discussed in the 
main text. 

We extract the static length scale from the above < 
Nc > data by using the function ^ Nc>, and we 

show in Fig. S5 the htting results of Eq. (2) and Eq. 
(4) of the main text. Similar to Fig. 3(b) and 3(c) of 
the main text, data collapse on the line defined by these 
two equations, indicating that the defined static length 
scale correlates well with the degree of supercooling and 
also the relaxation time. All the above results indicate 
that the definition of the densely packed particles has 
no qualitative influence on the results, and a reasonable 
definition is expected to obtain similar qualitative results. 


Results obtained from another definition of densely 
packed particles 


Discussions on some fitting parameters 


To study whether the definition of dense-packed par¬ 
ticles influences our results, we check our data by using 
another definition of them. As described in the main 
text, the statical distribution of the normalized Voronoi 
volumes is universal in the density range we studied for 
both LJ and WCA systems. Inspired by it, particles 
whose Voronoi volume satisfies {v— < v >)f<Ty < —1.2, 
are defined as the dense-packed ones. Based on this def- 


Next we discuss the values of A and extracted from 
fitting the T dependence of ^ with a function of Eq. (2) 
shown in the main text. We plot in Fig. S6 their depen¬ 
dence on density p for both LJ and WCA systems. Since 
the values quantitatively depend on the definition of the 
dense-packed particles, we obtain two sets of fitting pa¬ 
rameters, i.e., A\ (^ol) and A2 (^o2), which correspond 
to the definition of densely packed particles used in the 
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FIG. S 6: Density p dependence of the fitting parameter A 
(a) and (b) for both LJ and WCA systems. A1 and ^ol are 
obtained based on the definition of densely packed particles 
used in the main text while A2 and ^o2 are obtained based 
on the definition of densely packed particles used in the SM. 


FIG. S 5: (a) The rescaled length scale versus [ro/(r — 
Tq)]^. The black solid line is the relation of = [To/(r — 
To)]'^. (b) Relation between and [l/DZop(ra/ro)]'^. The 
symbols are the same as in (a). The black solid line is the re¬ 
lation of 5/^0 = [^/Dlogira/Ta)]^■ The definition of densely 
packed particles is shown in the SM text. 

main text and in the SM, respectively. In Ref. Q, the 
temperature T dependence of length scale (^') extracted 
from the structural order parameter can be described by 
= ^o\To/{T — , where d is the spatial dimension 

of the system. Note that the adjustable parameter 
used in their work is the order of the particle radius, but 
the fitting is poor for us to use this value as our hr 
Eq. (2) of the main text. As a result, it is meaningless 
to compare the value of A with 2/d reported in Ref. Q. 
We have checked that our ^g value roughly corresponds 
to the static length scale at the onset temperature (or 
a little below than the onset temperature) for each den¬ 
sity in both LJ and WCA systems. It has been reported 
that the temperature dependence of the static point-to- 
set length scale of the LJ and WCA systems collapse 
when the length is scaled by its onset tenmerature value 
and T is scaled by the onset temperature , although in 
a relatively small temperature range. In a larger range of 
supercooling this scaling collapse is found to break down 
in the hard spheres system [^. Our results may provide 
a more universal relation between the static length scale 
and temperature, as shown in Eq. (2) in the main text. 


It can be seen from Fig. S6 that the p dependence of 
these two sets of A and qualitatively similar. The 

value of A (both A1 and A2) decreases as increasing p 
for LJ system while this value has little dependence on p 
for WCA system. The difference of value A1 (or A2) for 
the two systems shrinks as increasing p and is about to 
disappear at p = 1.9. Similar behavior occurs in (or 
^o2). This is quite reminiscent of the different dynam¬ 
ics for LJ and WCA systems at different p. Our results 
indicate that the two fitting parameters A and are 
depending on the intrinsic properties of the system and 
are quantitatively impacted by the definition of densely 
packed particles. It should be noted that the range of the 
fitting parameters for different densities for both LJ and 
WCA systems is relatively small. However, we find that 
the fitting parameters are very sensitive to the fitting re¬ 
sults, and a little change of the value of fitting parameters 
can lead to obvious deviations from our data. Thus, the 
data collapse of Eq. (2) in the main text is not due to 
certain small values of A but due to the intrinsic system 
properties. 

Test of the generality of our method in other glass 
models 

To test whether it is general that the locally dense- 
packed structural environment is the structural signa¬ 
ture of slow dynamics, we study other two different mod- 
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els. The first one is the Wahnstrom (Wahn) system Q. 
It is a 50:50 binary Lennard-Jones mixture of spheres 
whose interaction parameters are very different from the 
Kob-Andersen model. The interaction parameters for the 
Wahn model are oabI^^aa = 0.916, ubb/ctaa = 0.833, 
^AB = CBS = CAA, and ms/rnA = 0.5. The potentials 
are cut and shifted at 2.5(Ta/3. In the following, simu¬ 
lation results will be described in terms of the reduced 
units ctaa, caa and mcr^^/cAA for length, tempera¬ 
ture, and time, respectively. We investigate the Wahn 
system with Nw = 3000 particles at the number den¬ 
sity pw = 1.297 in the temperature range T = 5.0 to 
T = 0.58. The mode-coupling temperature for this sys¬ 
tem is Tc Ri 0.56 [1,0. 

The second system we investigate is a 50:50 binary 
mixture of Nh = 3000 harmoni c sp heres (Harm) with 
the number density pu = 0.675 [lOj. All particles have 
the same mass m and the interaction potential is given 
by Vap{r) = 0.5e(l — rjuapY^ if ^a/3 < cFap and 0 other¬ 
wise. The interaction parameters are ctab/ctaa = 1.2, 
'^bbI^aa = 1-4. Length, energy and time are re¬ 
ported in units of ct^a, 10“^e and ^ respec¬ 
tively. The temperature range we studied for this system 
is T = 50 to T = 5.5 and the mode-coupling temperature 
is Te Ri 5.2 [l-liH. 




FIG. S 7: T dependence of the a relaxation time Tq of the 
Wahn system (a) and the Harm system (b). Solid lines are 
the fitting results by the VFT relation (see Eq. (3) in the 
main text). Here Tq is obtained from Fs{qp,Ta) = 1/e of the 
A part particles. The wave vector is equal to 7.7 and 6.28 
for the Wahn and Harm system, respectively. 

We also perform the molecular dynamics simulation 
with in-house GALAMOST software package [1| for the 
above two systems. Temperature is maintained by the 
Nose-Hoover thermostat and the integration time steps 
are 0.005 and 0.04 for the Warn system and the Harm 
system, respectively. Results are also averaged over eight 
independent runs for the two systems. 

We present the T dependence of Tq, and also the VFT 
fitting results for the above two systems in Fig. S7. 
Clearly, the VFT relation provides a good description 
of the results for both the Wahn and the Harm systems, 
especially at lower temperatures. 

To obtain the static length scales of these two systems, 
we also define 13% of particles with the largest local den- 




FIG. S 8: T dependence of the length scale ^ of the Wahn 
system (a) and the Harm system (b). Solid lines are the 
fitting results by Eq. (2) in the main text. 



FIG. S 9: Relation between and [l/Dlo(?(ra/ro)]'^ for 
the Wahn system (black squares) and the Harm system (red 
circles). The black solid line is the relation of = 

[1/Dlog{Ta/ra)Y^. Since the deviation from the VFT rela¬ 
tion (see Fig. S8) for the data points of the two highest tem¬ 
peratures of the Wahn system, we didn’t provide the data of 
T = 5.0 and T = 4.0 in this figure. 


sity as the locally dense-packed particles (this definition 
is the same as that in the main text). Then, we eval¬ 
uate the number averaged cluster size < Nc > formed 
by the densely packed particles. The static length scale 
is thus obtained from ^ We show in Fig. 

S8 the temperature dependence of the length scale ^ for 
both Wahn and Harm systems. Clearly, ^ increases as de¬ 
creasing temperature for both these two systems, which 
indicates that this kind of definition of ^ can also capture 
the static structure differences between high temperature 
and low temperature glass-forming liquids. On the other 
hand, the T dependence of ^ can be well fitted by Eq. (2) 
of the main text, indicating that the static length scale 
has a correspondence relationship with the degree of su¬ 
percooling (characterized by Tq/{T — Tq)) for both the 
Warn and the Harm systems. 

From Fig. S7 and Fig. S8, we can easily obtain the 
relationship between ^ and (see Eq. (4) in the main 
text). In Fig. S9, we plot against [l/Dlog{Ta/To)]^. 
Slight deviations from perfect data collapse can be seen. 
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which is because of the narrow range of the scaled ^ and 
the scaled , and little deviation of the data points from 
the fitting line in Fig. S7 or Fig. S8 can lead to obvi¬ 
ous deviations in Fig. S9. Nevertheless, the correlation 
between ^ and relaxation time is striking. Thus, we con¬ 
clude that the locally dense-packed structural environ¬ 
ment can also be an important structural signature of 
dynamics for both Wahn and Harm systems. 
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